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^ ' Abstract 



We investigate in some detail a "variational mass" expansion approach, gener- 



' alized from a similar construction developed in the Gross-Neveu model, to evaluate 

> ; 

V/^ \ the basic order parameters of the dynamical breaking of the SU{2)l x SU {2)r and 

fSj ■ SU{3)l X SU{3)r chiral symmetries in QCD. The method starts with a reorganiza- 

a^ : 

, tion of the ordinary perturbation theory with the addition of an arbitrary quark mass 

■ m. The new perturbative series can be summed to all orders thanks to renormal- 

I ization group properties, with specific boundary conditions, and advocated analytic 

continuation in m properties. In the approximation where the explicit breakdown 
pH I of the chiral symmetries due to small current quark masses is neglected, we derive 

ansatzes for the dynamical contribution to the "constituent" masses Mg of the u, d, s 
^ I quarks; the pion decay constant Ft^; and the quark condensate (qq) in terms of the 

basic QCD scale A-^g^. Those ansatzes are then optimized, in a sense to be specified, 
and also explicit symmetry breaking mass terms can be consistently introduced in 
the framework. The values of F.,^ and Mg obtained are roughly in agreement with 
what is expected from other non-perturbative methods. In contrast we obtain quite 
a small value of | {qq) \ within our approach. The possible interpretation of the latter 
results is briefly discussed. 
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1 Introduction 



Although QCD is certainly established as the correct theoretical framework underlying 
the strong interaction physics, a still challenging issue is to derive the low-energy proper- 
ties of the strongly interacting spectrum from QCD "first principle", due to our limited 
present skill with non-perturbative technics. At very low energy, where the (ordinary) per- 
turbation theory cannot be applied, Chiral Perturbation theory (ChPT) or extended 
Nambu-Jona-Lasinio models (ENJL) 0, |^ give a consistent framework in terms of basic 
parameters that have to be fixed from the data. Yet the link between those effective 
parameters and the basic QCD ones, the gauge coupling and the quark masses, remains 
largely unsolved. To have, typically, a determination of the pion decay constant and 
other similar low energy quantities directly in terms of the basic QCD coupling constant 
asi is clearly a desirable task. Lattice simulations are perhaps at present the only system- 
atic approach to such questions, and indeed provide, among many other things, a deter- 
mination of as [§] and also evidence for dynamical chiral symmetry breakdown (CSB) 0]. 
Yet a fully consistent treatment of the chiral symmetry on the lattice is still missing, and 
there are also inherent difficulties of dealing with truly dynamical "unquenched" quarks, 
light meson masses, and related problems. 

In the present paper, we investigate a new method to implement CSB directly from 
the basic QCD degrees of freedom, the quarks and gluons. More precisely we shall explore 
how far the basic QCD Lagrangian can provide, in a self-consistent way, at least the basic 
CSB parameters, namely the formation of dynamical quark masses, quark condensates, 
and the pion decay constant, in the limit of vanishing Lagrangian (current) quark masses. 
Such a qualitative picture of CSB can be made more quantitative by applying a new 
"variational mass" approach, recently developed within the framework of the anharmonic 
oscillator [^, and the Gross- Neveu (GN) model ||. 

Before developing our construction, let us remind that there are in fact two main 
phenomena which are believed to emerge from non-perturbative dynamics (and which are 
not fully understood from first principle), namely the Confinement and CSB respectively. 
Even if those are certainly intimately related in the full QCD dynamics, it is legitimate 
to consider, at least in a first approximation, those two issues separately. QCD is a very 
rich theory involving many different non-perturbative aspects, therefore any new method 
or model can most likely deal with a selected or simplified aspect of the full dynam- 
ics. After all, as far as the gross features of the chiral dynamics are concerned (that is, 
the existence of light pions and consequences in the form of low energy theorems 0), 
a picture without confinement but with CSB would not look very much different from 
the true QCD one. On the other hand, there is a large evidence that the confinement 
forces do play a more major role in the formation of heavier hadrons [jl^]. This relatively 
neat separation of the relevant scales is indeed one of the basic assumptions underlying 
several successful approaches to chiral dynamics, like the above-mentioned ChPT, ENJL 
models etc, which generally do not include the confinement properties in their framework. 
Somewhat closer to the present investigation, similar assumptions are often made in many 
dynamical CSB models, dealing with various levels of approximations in the treatment of 
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the Schwinger-Dyson equations for the Green functions relevant to chiral dynamics Q. In 
contrast, there exists a quite radically different attitude towards CSB in QCD, advocating 
that the responsible mechanism is most probably the non-perturbative effects due to the 
instanton vacuum or even more directly related to confinement [Q. However, even 
if the instanton picture of CSB is on general grounds well motivated, and many fruitful 
ideas have been developed in that context 0, as far as we are aware there is at present no 
sufficiently rigorous or compelling evidence for it. In any event, it is certainly of interest 
to investigate quantitatively the "non-instantonic" contribution to CSB, and we hope that 
our method is a more consistent step in that direction. 



Although the basic idea was already developed in refs 0, H (see also [|T5[ for a first 
acquaintance in QCD), we shall reformulate here the construction in some details for 
self-containedness reasons, and to take into account some crucial differences between 
QCD and the GN model. In addition we introduce explicit symmetry breaking mass 
terms consistently in the framework, which as we shall see play an important role in the 
determination of {qq). 

The method thus starts by considering an arbitrary quark mass term m added to the 
massless QCD Lagrangian, whose dependence is calculated in perturbation theory with 
an expansion parameter x, interpolating between the massive free theory, at x = 0, 
and the relevant interacting theory, at x = 1. This starting point is similar to the one 
developed a long time ago and implemented in various different forms in refs.||16|]-||27]|. 
There, it was advocated that the convergence of conventional perturbation theory may 
be improved by such a variational procedure in which the separation of the action into 
"free" and "interaction" parts is made to depend on some set of auxiliary parameters. 
The results obtained by expanding to finite order in this redefined perturbation series 
are optimal in regions of the space of auxiliary parameters where they are least sensitive 
to these parameters. Indeed, such regions are those expected to best approximate the 
exact answer, where there should simply be no dependence on the auxiliary parameters |. 
Moreover it recently appeared strong evidence that this optimized perturbation theory 
may lead to a rigorously convergent series of approximations even in strong coupling 
cases. In particular, the convergence of this variational-like procedure has been rigorously 
established in the case of zero and one dimensional field theories [0. 

An essential novelty however is that the construction in [0, ^ combines in a specific 
manner the renormalization group (RG) invariance of the theory together with analytic 
continuation in m properties: this, at least in a certain approximation to be motivated, 
allows to reach infinite order of perturbation in the parameter x, therefore presumably 
optimal, provided it converges. This leads to a set of non-perturbative ansatzs for the 
relevant CSB quantities, as functions of the variational mass parameter m and renormal- 
ization scheme (RS) dependence parameters as well, which can be studied for extrema 
and optimized. Quite essentially, our construction also provides a simple and consistent 
treatment of the renormalization, reconciling the variational approach with the inherent 



^ See e.g. ref. and the textbook by V. A. Miransky ||T^ for reviews and original references. 
^For a review on such issues and an extended coUection of original references, we refer to the textbook 
by E. V. Shuryak 0. 

^For variants of the optimized variational-perturbation theory ideas, not based on minimization, but 
leading in many field-theoretical applications to similar results, see also refs. p^, 
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infinities of quantum field theory and the RG properties. A comparison of the variational 
calculations with the known exact results for the mass gap in the 0{N) GN model [^] 
is very satisfactory in view of the fact that for arbitrary the (only known at present) 
second order non- logarithmic perturbative terms have been used in the optimization. For 
QCD, one might however argue that the truly non-perturbative effects are lying in the 
infrared domain, for which a complete non-perturbative treatment is clearly beyond the 
scope of the present framework. We shall see in fact that our RG invariant ansatz can 
be viewed as an (optimized) perturbation theory around a non-trivial fixed point of the 
RG evolution solutions, which in that sense avoids having to treat the infrared prob- 
lems inherent to non-perturbative QCD 0. The physical relevance of such a framework 
can be judged a posteriori, e.g. by comparing our estimates with the results from other 
non-perturbative methods (like the lattice simulations 0, or the spectral sum rules 
approach ||3^ typically) , but unfortunately there is at present no rigorous way of defining 
the intrinsic error and convergence properties of the method, at least in four dimension. 

The paper is organized as follows: in section 2, we explain and motivate in details 
the construction and derive the main ansatz for the dynamical CSB contribution to the 
"constituent" quark mass, in the exact chiral limit. A dynamical quark mass in QCD is 
not, strictly speaking, an order parameter, but it is interesting for at least two reasons: 
first it appears as a direct generalization of the rather intuitive concept of a (fermion) 
mass gap originating from dynamical CSB, as illustrated in many simpler models. Second, 
as we shall see, the concept and known properties of a gauge-invariant and infrared finite 
pole mass enters as a basic ingredient in the ansatz, and plays also some role in the 
determination of the genuine order parameters, (qq) and F^. Concisely, our construction 
will be shown to amount to: 

i) taking the pole mass, as expressed as a function of the Lagrangian (arbitrary) running 
mass parameter m{fi) in some appropriate renormalisation scheme (RS); 

ii) constructing an ansatz as a (contour) integral over the analytically continued mass 
parameter m{fi), which can be shown to formally resum the reorganized x series, with a 
specific contour and appropriate RS choices avoiding some singularities; 

Hi) finally optimizing the result with respect to m(/i) and the RS arbitrariness. 
In section 3 the construction is generalized to the case of a composite operator of arbitrary 
mass dimension n, and applied to the determination of the pion decay constant, F^. In 
section 4 we apply the construction to derive from similar arguments an expression for 
(qq) , introducing for this purpose explicit symmetry breaking mass terms. Summary and 
conclusions are given in section 5. We have also collected in Appendix A the relevant 
material on renormalization group properties, which fixes our normalizations and some 
definitions, and in Appendix B the analytic continuation technics used in the derivation 
of the main ansatzs. For completeness we also display in Appendix C the perturbative 
expressions which are used as basic ingredients in the different ansatzs. 



'*Quite recently, an interesting approach has been developed to determine the same CSB quantities, 
by using the exact "renormalization group flow" formulation [|34[ . Due to the central role played by RG 
properties, our present approach may be considered as having some similarities in aim with the latter 
one, albeit being totally different in the technics used. 
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2 Dynamical quark masses 



2.1 Preliminaries: QCD Lagrangian with variational mass 

Let us start with the bare QCD Lagrangian (for Uf quark fields qi of identical masses 

1 



QCD = —,G^:GI, + Y.qS^,D^ - mo)q^ (2.1) 
^ i=i 



where 



D% = d''q, + tgoTaG^,qr: (2.2) 

are respectively the gluon field strength tensor and the covariant derivative, a is the color 
index, i is fiavor index, Ta the SU{3)c generator, and mo, go the bare mass and gauge 
coupling constant, respectively. As is well known the QCD lagrangian (|2.1| ) possesses for 
mo = the additional chiral symmetry SUinf)^ x SUinf)^ x U{1)l x U{1)r, at the 
classical level 0. In what follows we only consider the SU{nf)L x SU{nf)ji part of the 
chiral symmetry, and for n/ = 2ornj = 3as physically relevant applications. 

Following the treatment of the GN model let us add and subtract a (bare) mass 
term, m-oZ^rii^i^j; to be treated as an interaction term, to the massless QCD Lagrangian. 
This is most conveniently done by introducing a new perturbation parameter x, according 
to 



Lqcd L^Qcoigo = 0, mo = 0) - mo ^ q^q, + LQcnidl xgl) + a; mo (2-3) 



- free 

vyo — '"-0 — ^) ~ '"-U 

i=l i=l 



where Lq^.^ designates the usual QCD interaction terms, i.e those proportional to go or 
9o (|2.1|) . This is formally equivalent to substitute everywhere in Lagrangian ( p.l|) 



mo mo (1 - a;); g^-^g^x, (2.4) 

and therefore in any perturbative (bare) quantity as well, calculated in terms of mo and 
go 0. Since the massless Lagrangian is recovered for x — > 1, mo is to be considered as an 
arbitrary mass parameter after substitution ( p.4| ). One expects to optimize physical quan- 
tities with respect to mo at different, possibly arbitrary orders of the expansion parameter 
X, eventually approaching a stable limit, i.e flattest with respect to mo, at sufficiently high 
order in x. 

One immediately encounters however a number of obstacles: first, and quite essentially, 
before accessing to any physical quantity of interest for such an optimization, the theory 
has to be renormalized, and there is an unevitable mismatch between the expansion in 



^We do not consider in particular the effects due to the breakdown of the axial U{1) symmetry 
associated with instanton phenomena [ p6[ . 

^ It is implicitly understood that the same coupling constant reparametrization, ( |2.4D , also affects 
e.g. the Fadeev-Popov terms necessary at the loop-level calculations, which we will not display explicitly. 



Note therefore that the substitution (2.4) does not spoil the gauge-invariance 
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X, as introduced above, and the ordinary perturbative expansion as dictated by the nec- 
essary counterterms, which at first spoils the whole picture. Moreover, for any physical 
quantity of interest in QCD one only knows at present the first few terms of the ordinary 
perturbative expansion in g'^, and the computational effort required to perform perturba- 
tive calculation in x is a priori similar. Finally, even if arbitrary higher order perturbation 
terms were available, it is easy to convince oneself that, at any finite order in the x expan- 
sion, one always recovers a trivial result in the limit m — (equivalently x ^ 1), which 
is the limit we are ultimately interested in to identify a mass gap (and more generally to 
identify spontaneous chiral symmetry breakdown, since an explicit Lagrangian mass term 
breaks the axial vector Ward-Takahashi identities explicitly |1T2|). 

All these problems can be circumvented at once, by advocating a certain resummation 
ansatz [§, exploiting the RG properties of the theory plus analytic continuation prop- 
erties in the arbitrary mass parameter m, which defines a candidate mass gap M{m) f\. 
With respect to ordinary perturbation theory this RG invariant ansatz is exact for the 
leading and next-to-leading dependence in m only, but the crucial point is that it resums 
as well the (reorganized) perturbation series in x. Moreover, taking the limit x — 1 in this 
resummed expression no longer gives a zero mass gap. The construction can be systemati- 
cally improved order by order (at least in principle, whereas in practice we obviously limit 
ourselves to the presently known highest perturbative orders for the relevant quantities). 
This is nevertheless hoped to be a sufficiently good starting point for further numerical 
optimization studies with respect to the (arbitrary) parameter m f\ For simplicity, we 
shall first illustrate the method by restricting ourselves in the next section to the first 
order of RG dependence. Generalization to higher orders will be discussed in the next 
sections. 



2.2 First order renormalization group mass ansatz 

To be meaningful our formalism should be applied to RG invariant and gauge-invariant 
quantities, in the QCD context. Therefore a reasonable starting point is to consider 
the pole quark mass, which is gauge-invariant |^ to arbitrary perturbative orders p^ . 



and infrared finite Using dimensional regularisation (with space-time dimension 



D = 4 — e), the one- loop expression for the quark mass, using the on-shell condition 
^ = mo, is easily obtained as 



mo [1 + 3Cpglj^^^^{mo)-^{l + -e + 0{e')) + 0{g',)] , (2.5) 

where Cp = {N^ - 1)/{2N,) for SU{N,) color (i.e. Cp = 4/3 for QCD), and r(x) the 
Euler Gamma-function. For e ^ 0, expression ( |2.5|) leads to the (renormalized) pole mass 



'' We will introduce several different definitions of a quark mass. To avoid confusion, we designate by 
m the Lagrangian "current" mass (or any mass trivially related to the latter) , while M (with possibly 
different indices) is used for dynamical mass candidates at different stages. 

* Such a situation is somewhat similar to many other variational approaches, where one often starts 
with a trial estimate, which one expects to be reasonably close to the exact result, before performing a 
variational calculation. 
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expression 



M[ ^ ^ (1 + 2^[| - ln(-)] + 0(a|) ), (2.6) 
no fi 



as expressed in terms of the MS mass m = m(/i) (/i = fxy/Ane ""^/"^^ 7^ ~ 0.577216.., 
where ^ is the arbitrary scale introduced by dimensional regularisation) , and coupling 



0:5 = g {fi)/{A7r). In (|2^) we used the MS mass counterterm Zm, defined as mo = 
and given from as: Z^{g^) = 1 - g^/{2'K^e) + 0(/). 



Our construction starts by first observing that RG properties gives us more information 
on higher orders than what appears in ( p^) , ( p. 61) : one-loop calculations actually provide 
the leading 1/e dependence of Z^ and of the counterterm for the coupling Zg to a// orders, 
as is well known: 

Z^^mo/m=(l + ^^r^, Z, = (^oV/K^=(l + ^^?')-^ (2.7) 

where 70 and 60 are the first order RG coefficient for the anomalous mass dimension and 
beta function respectively (see Appendix A for the normalization used). 

Now expression ( p.5|) using (|2.71 ) suggests the following form of a bare resummed mass 
expression (restricted at the moment to the first RG order): 

Ml = (2.8) 

[l-6or(e/2)(47r)^/2^2(Mi)-]^ 

in close analogy with the mass-gap formula of the GN model [§]. Eq. ( p.8|) is obviously RG 
invariant since expressed only in terms of tuq and g^. Introducing renormalized quantities 
m,g'^ from ( p.7|) , ( p.8|) is in addition finite to all orders, as easily checked, thanks to the 
recursivity in Mi, and reads explicitly 

TTl 

Ml = ^ . (2.9) 

[I + 2b,-gHnif)]< 

Eq. ( p.9| ) resums the leading log (LL) dependence in rh of the usual pole mass but, ev- 
idently, does not give its correct next-to-leading dependence in m (nor even any of the 
"purely perturbative" (non-logarithmic) finite corrections, such as the 2/3 factor appear- 
ing already at first order in ( ^.5[ )-( P^ ). Actually, the renormalized form ( |2.9| ) can 
alternatively be rigorously derived as a particular boundary condition imposed on the 
usual RG running mass. Indeed, consider the general solution for the RG evolution of the 
running mass, 

m(/i') = m{fi) exp{- f dg^^wr^} (2-10) 

f3{g) 

in terms of the effective coupling g{fi), whose RG evolution is given as ^dg{fi)/dfi = [3{g). 
Solving ( |2.10| ) with the following "fixed point" boundary condition: 



M = m{M), (2.11) 

one easily obtains (to first RG order) 

Ml = ^ ^ , (2.12) 

[l + 26o^?2(/i)ln(Mi)]^ 
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which is nothing but eq. ( p.9| ) for fi = fi. In particular, the recursive form of ( p.l2| ) with 
respect to Mi can only be obtained thanks to the specific boundary condition in (2.11). 



Note that, although { \^.9\ ) can thus be obtained solely from the condition ( |2.11| ) imposed 



on the "current" mass m{fi), it has the usual properties of a pole mass ^ This coincidence 
between the pole mass M and the current mass m{M), evaluated at the scale = M, 
is of course only an artefact of our crude approximation, neglecting at the moment the 
non-logarithmic perturbative corrections. Although the first order expression (|2.9| ) having 
those properties could be guessed directly, the advantage of the latter derivation is that it 
gives a sound basis and precise link with the renormalization group behaviour. Indeed, as 
will be derived in the next section, one can still obtain an exact solution of ( p.lO|) with the 



boundary condition (|2.11| ) at the next order explicitly, which is by construction entirely 



consistent with the next-to-leading RG dependence, to all orders. We shall also include 
later on the non- logarithmic perturbative corrections, necessary to make contact with the 
genuine pole mass and to define a more realistic "mass gap" ansatz. 



Now, the most important property of the simple-minded expression (|2.9|) , which is 
also shared by the more realistic ansatzs to be later derived, is that it is non-zero in the 
chiral limit, m — > 0. To see it explicitly, first rewrite identically (p.9|) as 



Mi(ln(Mi/Ajg5))^ =m (2.13) 



where for convenience we simply introduced the RG invariant scale J^jjg = /ie ^''of (at 
first RG order), and the scale-invariant mass m = fh{2bog'^) ^^o . For fixed g"^, the chiral 



limit is m ^ in (|2.13| ). Now, (|2.13|) may be seen as a function m(Mi), and requiring its 



inverse Mi(m) to be defined on the whole physical domain < m < oo, and to match the 
ordinary perturbative asymptotic behavior for m —>■ oo, implies Mi{m 0) —>■ A;^0. 
This property of ( p.9|) is in contrast with the "one-loop improved" expression of the mass 
(obtained by replacing Mi — > m in the denominator of (|2.9|) ), which would only give Mi — *■ 
for m ^ 0. In other words, while ( p.8|) -( p^ are perturbatively compatible with the 
usual (one- loop) RG behaviour of the current mass, a specific choice of boundary condition 
for the RG running has selected a non-trivial "vacuum", by changing the behaviour of 
the whole perturbative series in the chiral limit fh ^ 0. Alternatively, one may derive the 
result Ml = A.j^ in the chiral limit as follows. From (|2.13|) , defining 

F = H-^].M^]-2L,^F, (2.14) 

on the physical branch, < m < oo, the latter relation can be inverted as 

e^F^ = . (2.15) 



^In particular, Mi in ( ^.12| ) is scale invariant, in contrast with m(/i): i.e. an arbitrary change in 
in the current mass m{^) is compensated, to all orders, by the In/x dependence in the denominator of 
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(2.12). (2.12) is also clearly gauge-invariant, as the pole mass should be. 

There is a priori another possible solution of ( ^.15 ) in the chiral limit: Mi — > when m ^ 0. 



However, it is easily seen that Mi{rh) has branch points at Mi — Ajj-g and Mi = e Ajj-g. The 

branch giving the trivial solution is only defined for < |m| < (7o/25o)''''/2boe-To/2boyY— < A^^g^, and 



is therefore not compatible with the asymptotic perturbative behavior of (2.9) for m oo 
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Now it is a simple algebraic exercise to see that ( |2.15| ) gives F{m) as an expansion in 
powers of (m/ Ajjg)^'^'''^^"'°\ for small m: 

F{m 0) ^ (-— ) -0 [1 - ^(-— ) -0 + . . . ], (2.16) 

which together with (|2.13| ), i.e. Mi = m F~^'^°^'^^°\ immediately gives Mi foi' 
m — > 0. While the latter derivation may be unecessarily sophisticated at first order, it 
has the advantage of being the one easily generalizable to higher order, where we shall 
also find a non-zero mass in the chiral limit, M = const . Ajjg, but with a no longer trivial 
proportionnality constant. More importantly, it is also the most convenient procedure to 
analyse the chiral limit once having established an expression analoguous to ( |2.9| ), but 
obtained directly from the variational perturbative (bare) expansion (as defined by the 
substitution ( |2.4| )), which we address next. Our aim is to obtain a variational "mass gap" 
where the non-trivial chiral limit property of ( |2.9[ ) is preserved, while at the same time 
providing us with a systematically improvable ansatz, in accordance with the "variation- 
ally improved perturbation" principle. 



Let us thus proceed with the x parameter expansion, performing in (|2.8| ) the substi- 
tution (|2.4| ). Unfortunately, as announced in the last section, it is easily checked that the 
resulting expression is no longer finite: after the introduction of the usual counterterms 



in (|2.7| ), some of the divergent terms in the denominator of eq. ( p.8|) no longer cancel, 
except for a; — >^ 1, in which case one only obtains the trivial result Mi —>■ 0. This is 
not much suprising since the precise all-order cancellations of divergences in eq. ( p.8|) , 
resulting from the above RG properties, have no a priori reasons to be compatible with 
the peculiar modification of the basic Lagrangian mass and coupling as implied by the 
substitution ( p.4|) . The way out is to resum the x-generated series, denoted by Mi(a;), in 
a different manner [Q: by analytical continuation in x one can find an adequate integra- 
tion contour, resumming exactly the series Mi(a;) in the a; — >■ 1 limit. This is explained 
in full details in Appendix B. The net result gives an ansatz (still at the first RG order 
approximation) : 

1 f dv V mo 



where v = q{l — x) has been introduced as an appropriate change of variable to analyse the 
X ^ 1, g — s> oo limit of the g^'^-order expansion of Mi(x), the contour is counterclockwise 
around the negative real axis, and we have defined 

Mv) ^l-bo gl r[-] [movY^ ifoY^ (2.18) 



whose expression is directly dictated from the denominator of eq. ( [^.8] ) (see Appendix B 
for details). The crucial point is that ( |2.17| ) is now finite to all orders, while giving a 
non-zero result. More precisely, after renormalization eq. (|2.17D becomes 



1 r rrj 

Ml = i dve" ^ (2.19) 

2«7r J f{v)^''o 
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where the renormahzed function f{v), related to fo{v) above as Q 



/o(l) = ZJil), (2.20) 
satisfies the (finite) recursion formula 



_2 , ,m V „_au- 



/ = l + 26orin(^/ ^). (2.21) 



Introducing a last convenient change of variable, 



eq. (|2.19| ) takes the form: 



where 



/ 1 

^^2b^ (^II) 2bog\M,) ' ^^-^^^ 



M, = ^idyey/-^'^, (2.23) 



m" = (m/A^) (260^^)"^"/^'"' ; y = m"v (2.24) 



are just the conveniently rescaled, dimensionless (scale-invariant) "mass" parameter and 
integration variable, respectively. After those different manipulations, it is easy to see 
that one can use (|2.14|) to express F in (|2.23|) as an expansion for small y, using (|2.16|) 
with the substitution {rh/Ajjg) —>■ y: 

Fiy ^ 0) ~ ?/ ^0 (1 "iy 70 + Oiy ) ). (2.25) 

7o 

This is nothing but a manifestation of the dimensional transmutation mechanism in the 
m — > limit of QCD, although in a rather unconventional form: the function F, defined 
in ( |2.22| ) as a function of the running coupling becomes, for m — 0, a function of the 
(rescaled) mass only, y = m"v. The original expression depending on {(7^, ffi} now depends 
on {Kjjg, m"} (where m" is arbitrary), and we are ultimately interested in the chiral limit 
m" 0. From ( |2.25D it is clear that ( ^.23| ) has a simple pole at y ^ 0, with residue 
giving the announced non-trivial solution. Mi = Ajjg. 

The previous construction therefore shows that summing the x series for the varia- 
tional (bare) expansion and renormalizing, obtaining the contour integral in ( p.l9| ), is 
equivalent to performing the following steps: 

-i) take the renormalized RG solution of (|2.10| ), M(m), with the condition ( |2.11|) ; 



-ii) performing in M{rfi) the substitution fh — >■ ffiv, and integrating the resulting expres- 
sion around the cut negative real axis, with a specific weight §{dv/v) e"" as in eq. ( p.l9|) Q- 
Having this well-defined integral transform we can work directly with renormalized quan- 
tities and define formally an ansatz at a given order as an appropriate generalization of the 

The simple relation in ( |2.20|) , involving the counterterm Zg only, is an accident of the first RG order 
approximation. At higher orders the relation is much more involved |q|. 



Actually at the next order, solving ( 2.10 ) and taking the integral afterwards is not strictly equivalent 
to the derivation starting from the bare expressions. However one can still show |^ that the two derivations 
are simply related by a particular renormalization scheme change. 
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integral in ( |2.19| ), provided however that one does not encounter extra singularities in v. 
For simplicity let us postpone the question of the possible existence of extra singularities 
until section 2.4, where this potential problem is addressed in details. 



It is now straightforward to introduce the necessary non-logarithmic perturbative cor- 
rections to the (purely RG) above results: this is consistently done as 

M[{m") = ^ ldyeyl-"^{l + (|)-^ + 0(i,) ) (2.26) 



without changing anything in the contour integral properties, except that (|2.26|) no longer 



has a simple pole behavior at y —>■ 0. Eq. ( |2.26|) is nothing but a specific integral over the 
pole mass: namely taking the integrant for v = 1 and expanding for small by using 
relations (|2.21|) and ( p.22|) , one explicitly recovers the usual pole mass expression with 
first order non-logarithmic perturbative correction. 

One may now optimize ( |2.26| ) with respect to m", by performing (numerically) the y 
integral. In fact, numerical integration is not mandatory: since what we are interested 
in is the behavior for m" ^ 0, it is equivalent to look at the properties of eq. ( |2.26| ) for 
y ^ 0. An expansion of ( p.26[ ) near the origin is provided from Hankel's formula. 



1 /■ , " (m V^" 

— I dyey'^ y- = , (2.27) 

ziTT J i [—a\ 

where the different powers a resulting from the expansion at arbitrary order near the 
origin have the form a = —1 + P^, with p = — 1, 0, ■ ■ ■ integer 

We stress that m" (equivalently, m) is meant to be an arbitrary but implicitly small 
parameter: even if we ultimately seek for optimal values ^rn'^p^ as best approximations to 
the limit m" 0, at intermediate steps m" ^ (i.e m 7^ 0) implies that the axial vector 
current Ward-Takahashi identities are explicitly broken. Eventually a small, physically 
acceptable explicit breaking (like the PCAC hypothesis 0) may be considered, where m^p^ 
ip^o-pt) could be tentatively interpreted as the actual physical quark masses. But that is 
largely in contradiction with the basic principle adopted here, according to which m" 
itself has no physical meaning, whereas only the optimal value of {m") in (|2.26|) has. 
Indeed, in the (ideal) situation where the optima would be really fiat, m" would obviously 
not be well determined. It is therefore much more preferable to find a mean, eventually 
an approximate one, to reach the exact limit m" 0. The chiral Ward identities are then 
recovered and we can define the mass gap strictly in the chiral symmetric limit. We shall 
see in section 4 how to implement consistently an explicit symmetry breaking, physical 
mass term, independent of the arbitrary variational mass parameter m". 



Note that (2.26) with (2.27) has some similarities with a Laplace-Borel transform: indeed the 
analytic continuation in m is almost equivalent (for m — > 0) to a continuation in the coupling g. The 
precise connection, and the link with the renormalon singularities | ]39| associated with the usual Borel 
transform will be discussed in more details elsewhere 40|. We shall however come back in the next 
sections on the structure of singularities of the integral ( 2.26| ) encountered in our framework. 
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2.3 Second order dynamical mass ansatz 

At the next RG order, the previous quahtative picture remains essentially unchanged, 
except that the derivation is somewhat more involved. The solution of (|2lO| ) with (|2lTl) , 
using the two-loop RG coefficients ffl, given in Appendix A, is 



Mo 



where / = g^/g'^{M2) satisfies 



m f ^''0 



(2.28) 



/ = l + 26orin— + 7^rin 

/i &o 



/ 



(2.29) 



(note in (|2.28|) - (|2.29|) the recursivity in both / and M2), and the non-logarithmic pertur- 
bative corrections are easily included as 



MJ, 



M2 (1 + -70^ + 



2 r K -g' 

^(4712)2/2 



+ 0{g' 



(2.30) 



where the complicated two-loop coefficient K was calculated exactly in ref. |^3| and is 
given explicitly in Appendix C. Eq. (|2.30|) defines the pole mass including two-loop 
non-logarithmic corrections, and can be easily shown in addition to resum the leading 
and next-to-leading logarithmic dependence in m to all orders (see Appendix A). The 
contour integral generalization of ( p.26| ) is obtained, as explained before, after making the 
substitution m — > mv in ( p.2^ )- (|2.3(]| ), as 



(m",a) 



A' 



MS 



2in 



dy 



F{yY[C + F{y)] 



B 



F{y) F{yy 



(2.31) 



where ^jjg is now the RG invariant basic scale at two-loop order, in MS 



IIP 2bo9- 



'1 U il 



(2.32) 



^0 



and the dimensionless (scale-invariant) arbitrary mass parameter reads 

1 B 



m" = 



, m 
• A 



2^ [2h^g 



2l-7o/(26o) 



bo 



In (|2.31|) F now satisfies the recursive relation 

Fiy) ^ \n[y] - A ln[F(y)] - (B - C) \n[C + F{y)], 
where A, B, C are expressed in term of RG coefficients as 



^ ^ 7o 7i 
26o 26i ' 



(2.33) 



(2.34) 



(2.35) 



^*In ( 2.29 ) the recursivity operates on Af2, which by definition o nly d epends on the RG evolution, 
therefore to be distinguished from the (pole) mass expression in (2.3C). 

^^With our prescritions one should keep in mind that ^jjg depends on ny, as usual p5[ |. 
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explicitly given in Appendix A. 



In ( |2.31D we have also introduced an extra parameter a, taking into account changes 
in the arbitrary renormalization scale, according to fl —>■ a fi. Indeed, even at the first RG 
order there are infinitely many ways of introducing the non-logarithmic perturbative cor- 
rections to ( p.23| ): for instance, we might have introduced the 0{1/F) terms via another 
definition of F: F' = F{1 + 0{1/ F)) in place of (|2.14|) , up to second order terms. This 
non-uniqueness is basically what is parameterized by a, and accordingly the perturbative 
coefficients M.i in ( p.31|) have a logarithmic dependence in a, simply dictated order by or- 



der from the requirement that ( p.31|) only differs from the original MS scheme expression 
by higher order terms: 

A<i(«) = ^(^-ln a), (2.36) 

•^'^''^ = (2^^ (4^ ^ ^ ^""^ ~ ~ f + ■ ^^-^^^ 

Since we shall however continue expression ( |2.31| ) to the (non-perturbative) region of 
m" — >■ 0, in the manner described in the next sections, the a-dependence will eventually 
exhibits non-trivial extrema, and it is thus sensible to optimize the result with respect to 
a Q, since the unknown exact result would not depend on the arbitrary renormalization 
scale /i. Note that a ^ 1 corresponds to a change of renormalisation scale only, while keep- 
ing the renormalization scheme fixed, i.e. g"^ = g'^ifi) (equivalently Ajj-g) and rh = m{fi) 
fixed (to MS). Actually the situation is slightly more complicated, since at second per- 
turbative order there are other possible changes of renormalization prescriptions than a 



simple change of scale, which do affect expression ( p.31|) and should therefore be taken 
into account in principle. This turns out to be an important aspect of our analysis and 
will be addressed in details in the next section. 



For small y, F{y) from (|2.34|) has the expansion: 



F{y ^0) A{u- + 0{u^)); u = A-^C^y^, (2.38) 



again implying that ( |2.31| ) would give a simple pole at y — > for vanishing perturbative 
correction terms, A4i = = 0, with residue {2C)~^ . Note that all the results of section 
2.2 are consistently recovered by taking 6i = 71 = and neglecting the non-logarithmic 
perturbative corrections of 0{g^) in the different expressions above, ( p.28| )- (|2.37| ). 



It is still possible to generalize to the next RG order the previous derivation, since in 



QCD the three- loop coefficients 62? I2 are known ^g], although the construction becomes 
quite cumbersome. In fact such a generalization is not worth doing, since one can in fact 
choose a renormalization scheme in which 62 and 72 are set to zero, as well as all higher 
order coefficients. This is examined in more details in the next section. 



This procedure indeed gave very good results in the GN model, where in particular for low values 
of N the optimal values aopt are quite different from 1. 
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2.4 Renormalization scheme changes and generahzed ansatz 

As it turns out, the simple picture emerging from the last section is unfortunately ques- 
tionable as one realizes that, strictly speaking, formula ( p.31| ) has extra singularities in 
the y plane, in addition to the cut on the negative real y axis, implicit in its derivation 
and in ( p.27| ). After all, this specific contour was suggested by the known properties of the 
GN model, and it is not surprising if the analytic continuation in m properties are more 
complicated in QCD. To begin, as is clear from (|2.31|) there is an extra branch point at 



F = —C = — 6i/(26q), which was not present in the first order ansatz ( p.26|) , and simply 
corresponds to the first non-trivial fixed point of (3{g) located at = — 6o/&i- Fortunately 
it is harmless in the present QCD case, since 61 > implies that it is located along the 
already cutted F < real axis 0. 

Less trivially, the zeros of dy/dF also give from ( p.34| ) extra branch cuts in the y plane, 
starting at the (complex conjugate) points 

Vcutili) = e^(^^)F^(7i)[C + F(7i)](^-^) ; (2.39) 

where 

^(7i) = ^(-:^)[l±(l-4^)^/T (2.40) 

However, as can be seen the actual position of those branch points do depend on the 
scheme via the second coefficient of the anomalous mass dimension, 71 (the precise RS 
dependence of the latter is given in eqs. ( [A .61 ), ( [A .81 ) of Appendix A). In the original MS- 
scheme, F{'yf^^) gives extra cuts starting at Re[ycut\ — 0.34 (0.24), Jm[?/c«t] — ± 0.74 
(±0.76) for Uf = 2 (jif = 3) respectively. Note that when looking at the limit of interest, 
y — 0, using the expansion in ( |2.38D , one never sees those extra cuts: at any finite order 
the series expansion only has the cut on the negative real axis according to eqs. ( |2.27| ) and 
( ^.381 ). These singularities may after all be an artifact of our extrapolation to very small 



m" of a "perturbative" (although resummed) relation, eq. ( p.34|) . Anyhow, the point is 
that using the expansion near the origin is invalidated if there are extra singularities lying 
in the way with Re[ycut] > 0, since it would lead to an ambiguity of 0{exp{Re[ycut]/^")) 
for m" — * 0, in the determination of the integral (|2.31|) 0. 

The way out is thus clear: if there exists values of 71 which move those extra cuts 
away (or which are such that they start at Re[ycut] = 0), the expansion around the origin 
in y is legitimate. Defining 

7^ = 71 + A71 , (2.41) 

^^In the GN model case the similar cut was more troublesome, since bf^ < 0, and we used a specific 
Pade approximant constructio n [p| avoiding the cut. In a different context, this non-trivial fixed point 
has also been studied recently [l47| in connections with renormalon properties. 

^^This is quite similar to the renormalon ambiguities |p9|. Although the usual renormalons are not 



explicitly seen here due to the appearance of a mass gap (i.e., within the mass gap ansatz (2.31), by 
construction there is no integration over the Landau pole region), the way in which the singularities in 
y appear, namely in a resummed expression relating a "reference" scale Mdyn — ^jjs infrared 
scale m" ~ 0, may be viewed as reminiscent from the renormalons. An essential difference is that in 
the present case those singularities occur in the analytic continuation of the mass parameter rather than 
the coupling, and that it is possible to move those singularities away by appropriate RS change, as is 
discussed below. 
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one easily finds that Re[ycut\ — for A71 ~ 0.00267 (0.00437) for Uf = 3 (nf = 2) respec- 
tively Therefore we can adjust a correct 7^ by performing (perturbative) RS changes. 
As explicit from eqs. ( [A. 61 ), ( |A.8D of Appendix A, this can be done either by a first order 
change m. g g or m m' (or both). One thus considers the most general RS change 
(but restricted to the second perturbative order, which is sufficient for our purpose) and 
optimize with respect to this new arbitrariness. (As a side remark, we mention that one 
would also have Re[ycut] < with 71 = 7^^'^, for > 5). Note that the removing of 
the unwanted singularities is only possible at the two-loop RG order ansatz, eq. (|2.31|) , 
due to the first occurence of RS arbitrariness at two- loop order via 71, in MS schemes. 
More precisely, for the first order ansatz, eq. ( p.26| ), an extra singularity occurs now with 



Re[ycut\ = -Re[exp [—70/260] (—70/260)'^°^^'"'] > 0, which accordingly cannot be removed 
by a RS change. 

We obtain after simple algebra the generalized dynamical mass ansatz in the new 
(primed) scheme, in terms of the arbitrary RS change parameters Ai, A2, -Bi, B2 defined 
in Appendix A): 

^MS " 2- f \F'iy))^'[C + F'iyr^'+ F'{y) + F^(y) ^' ^'-"'^ 

with 



2h 



'0 



M'2{a,A,,B,,B,)= M2{a)-{A^ + B,)^^-'-^-^ (2.43) 

ZOo 40n 



M\{a) (52-7o5i) 



and 



^MS = ^^P\^]^-MS (2-44) 

(where relation ( p.44| ) is exact to all orders |^) and F' and m" have specific expansions 
in terms of the original MS quantities F and m" that we will not need explicitly here. 
We also impose a further RS choice, 

6; = 0; 7^ = 0, (2.45) 



which, according to relations (|A.7| )- (|A.8| ) in Appendix A, fixes A2 and B2 uniquely in 



terms of Bi and A71 and guarantees that the definition of A;^ in (|2.32| ) is unaffected. 



apart obviously from ( p.44|) . In what follows we express all results in terms of the original 



^Nis scale. 

2.5 Pade approximants and numerical results 

It is worth emphasizing at this point that our purpose is not to find a particular RS 
choice, which would "best fit" the expected order-of-magnitude result for the (dynamical) 

Actually Re\ycut\ is a (semi)-periodic function of A71. We only consider solutions of i?e[j/cut(A7i)] ~ 
nearest to the original MS value (A71 = 0) for nf = 2, 3. 



Although e.g. A2 does not appear explicitly in (2.42), we nevertheless need a prescrition to fix it, 
since it appears at the second order general RS change. 
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constituent mass (or similarly for and m{qq) considered in the next sections). Indeed, 
it is very likely that with so much RS freedom at disposal, one could make eq. ( |2.42| ) 
fit almost whatever values one wishes. In contrast, what we are seeking is the flattest 
region in the arbitrary RS parameter space, in the "Principle of Minimal Sensitivity" 
(PMS) sense. One soon realizes however that our extension of the PMS defines a 
rather complicated optimization problem: one has in principle to find the flattest possible 
extrema of ( p.42|) , in the three independent parameter space {m", a, Si}, where in addition 
A71 is constrained to give good analyticity behavior of ( p.42|) . 



Fortunately, one can study this problem within some approximations, which we be- 
lieve are legitimate. As above explained the ansatz in (|2.31|) (or ( p.42D ) is already optimal 



with respect to m at m" = 0, by construction, for vanishing pertubative non-logarithmic 
corrections Aii = 0, and in this case the optimal result is the simple pole residue. Due to 
the non-logarithmic, purely perturbative corrections (which are at present only known to 
second order for the mass), this simple picture is lost, but accordingly one may assume 
that the resulting expansion for small m" is "as close as possible" to an optimum (as 
would be the case if the sequence of approximations obtained by considering increasing 
orders of the x expansion could be proved to converge, like in one- dimensional mod- 
els p^). Accordingly, we will define the m" limit of ( |2.42| ) by a relatively crude 
but standard approximation of those perturbative corrections, rather than performing a 
numerical optimization with respect to m". Indeed, as discussed at the end of section 
2.2, it is in addition physically motivated to reach m" —>■ 0, since the axial vector current 
Ward-Takahashi identities are recovered and a non-zero result signals the spontaneous 
chiral symmetry breakdown. 

The approximation we are looking for is certainly not unique: given the ansatz (|2.42|) , 
one may construct different approximant forms leading to a finite limit for m" — *• 0. 
We shall demonstrate the feasibility of our program in the simplest realization: since the 
resulting expression will anyhow be optimized with respect to the RS dependence (entering 
any such approximants via the RS dependence of we assume that it largely takes 
into account this non-uniqueness due to higher order uncertainties, in the standard PMS 
sense. The latter assumption is also supported by the results in the GN model [§, where 
several different approximants were tried and compared with the known exact results. 



Pade approximants are generally known to greatly improve perturbative results 
and in most cases have the effect of smoothing the RS dependence. We shall consider the 
following Pade approximant for the purely perturbative part of the integrant in ( |2.42| ): 



PadeiF)^^^±^^ - M, (2.46) 
^ ' F + p{a,Bi) F^O p(a)' ^ ' 

which by construction restitutes a simple pole for F — in ( p.42| ). Matching its pertur- 
bative expansion for F 00 to the one in ( |2.42| ) one obtains 



M^^'^'ia, A71, Bl) = Aj^ i2Cr^ a exp{^} 



A^2(a,A7i,5i) 



(2.47) 



'^■^ For an extended discussion of the PMS motivations, we refer to the original Stevenson's paper, 
ref. ra. 
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When defining ( p.46| )-( p.47| ) one should be careful not to introduce new poles in the F > 



{y > 0) plane. This simply implies that there are some constraints on the possible values 
of (a, Bi) which, if we are lucky, do not invalidate their optimal values obtained from 
minimization. 

We have performed a rather systematic study of the possible extrema of the Pade 
approximant formula ( p.47|) for arbitrary a, Bi, with A71 fixed such that the extra cuts 



start at Re[y] ~ 0. We did find an optimal region with respect to a, in the sense that it 
minimizes the second derivative at the maximum 0. The region of the parameter space 
in the vicinity of that extremum is illustrated in Fig. 1, as function of a for different 
values of Bi. The plateau region has been determined more accurately using a numerical 
steepest descent method, requiring a minimal curvature. (For instance for Uf = 2 it 
corresponds to a^pt — 1.42, Biopt — 0.12). Our results are also summarized in Table 1 
(see section 5). The optimal values of the parameters a, Bi, A71 are consistent with the 
further requirement that ( |2.46| ) has no poles at F > 0, as it should. Explicitly we obtain 

Mi^'^'iopt) ^ 2.97 A^(2) (2.85 A^(3) ) , (2.48) 

for = 2 {uf = 3) respectively. 

We shall see in the next section that we can relate the pion decay constant F^, to Ajjg{2) 
or Ajg^(3), from which one may eliminate to obtain an evaluation of M^'"'^'^ in the 
chiral limit. Anticipating on these determination of A;^, eq. ( |2.48| ) with ~ 92 GeV 
gives (still in the pure chiral limit) 

M^'"^'^{opt) ~ 500 (447) Mel/ (2.49) 

for Uf = 2 [uf = 3). 

Finally, for a useful comparison, we also give the value of (|2.47|) in the MS scheme (which 



may be thus considered as optimized with respect to m" but in a fixed RS): 

Mi^'^'iMS) ~ 0.99 Aj^i2) (1.00 Aj^i3) ) . (2.50) 



Note however that according to the previous discussion, the result in ( |2.50| ) are a priori 
plagued by an ambiguity associated to the bad-placed extra singularities in the origi- 
nal MS scheme. This may however give a qualitative idea of the effects of the above 
optimization with respect to the RS choice. 

3 Generalization to Composite operators: F-j^/A 



MS 



3.1 RG-invariant ansatz for a composite operator 

We shall now derive an ansatz similar to eqs. (|2.31| ), (|2.42|) for the pion decay constant 
F^. The main idea is to do perturbation theory around the same RG evolution solution 



^^Note that the second derivative at the extrema points is a good quantitative estimate of the flatness 
criteria even in a multi-dimensional parameter space: the intrinsic curvature at an extremum of the 



hypersurface as deflned from e.g. the approximant ( ^.47 ), is proportional to the product of its second 
derivatives with respect to the different parameters. 
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with the non-trivial fixed point, as specified by the function F in ( p.34| ), with perturbative 
corrections specific to obviously. One should first identify from a relation where 
the previous construction can be best generalized. A formal definition which suits all 
our purposes is the well-known low-energy expansion of the axial vector-axial vector two 
point correlation function |5^ 



I J d'qe''^-^ < 0|T 4(x)A^(0)|0 >^ 5'^g,,Fl + 0{p,p,) . (3.1) 

In (|3.1j ) T is the time-ordered product and the axial quark current, = (97^75 A*g)/2, 
where the A*'s are Gell-Mann S'?7(3) matrices or Pauli matrices for nj = 3, Uf = 2 
respectively. The non- vanishing of expression (|3.1|) implies CSB: in other words F^^ is 



to be considered as an order parameter [51, E2, B3|. Since F^^ is expressed in terms of 



a gauge-invariant and RG invariant composite operator in (|3.1| ), one can apply a rather 
straightforward generalization of our ansatz exploiting in particular the RG resummation 
properties. The perturbative expression of (p.l|) for m 7^ is known to the three-loop 



order and is given explicitly in MS scheme in Appendix C. 

With a little bit of insight, a generalization of ( ^.28[ )-( p.30[ ) to a composite operator 



of naive mass dimension n (depending only on g and m), can be written as 

^ ^^j^^a^^cTfy^ ^ f f^ ^''^' ^ ^ 



in terms of F as defined by eq. ( 2.34| ) (for y/m" = v = 1), and where 6, a, (3 are fixed by 



matching the perturbative expansion in a way to be specified next. In practice, we shall 
only consider (|3.2| ) for the relevant cases of n = 2 and n = 4, corresponding to F^ and 
m{qq) respectively. 

Formula ( p.2|) as it stands is not yet our final ansatz and necessitates some comments. 
Apart from the trivial powers of n dictated by dimensional analysis, a rather obvious 
difference with the mass formula is that the composite operator being not a Lagrangian 
term, its perturbative expression starts at the one-loop, but zero*^ order in g'^. To acco- 
modate this fact with the correct RG properties, the expansion of (3^) starts with a l/g"^, 
due to the extra 2boF factor. The l/g"^ first order term anyhow cancels after a necessary 
subtraction which we discuss now |^. 



A more essential difference with the mass expression is that ( |3.1| ) is not finite even 
after mass and coupling constant renormalization. One can define a new, finite quantity 
after an adequate subtraction. In the present case, since our construction starts from 
perturbative expressions, this turns out to be nothing but the usual way of renormalizing 
a composite operator. Accordingly the prescrition is (perturbatively) well-defined [^] 
and unambiguous for a given RS choice. What is lost, however, is a part of the predictive 
power: unlike the mass case, a consistent treatment of the subtracted terms (i.e. respect- 
ing RG invariance) implies that the unambiguous determination of the 1 / perturbative 
terms in the final ansatz necessitates information on the {n + 1) order of perturbation 
theory. 



^^Note that the form and properties of (3^) can alternatively be rigorously derived by following a 
construction starting from bare quantities, similar to the one explained in for the vacuum energy of 
the GN model. 
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In the two cases here considered, and rh{qq) in the next section, this renormahzation 
actually reduces to a simple subtraction of the operator [m" 1] This implies, from RG 
invariance, residual finite subtraction terms, affecting the perturbative correction terms 
a and /3 in (^.21). Accordingly (|3.2|) is modifed to 



where it is easily shown, using RG properties, that the finite subtraction function H{g 
X^i^o HiQ^^ is determined perturbatively order by order (in the MS scheme) from 

[nim{g) + - m-^] H{f) = f • (3.4) 

In ( |3.4|) Ci((7^) is given by the residue of the 1/e term in the perturbative series expansion 
for the relevant quantities, as given in (|C.3|) and ( |C.10| ) respectively for F| and m{qq). The 



consistency of our formalism is checked by noting that the expansion of ( |3.3| ) in powers of 
g"^ do reproduce correctly the LL and NLL dependence in fh of the perturbative expansion 
of the composite operator to all orders, as well as the perturbative non-logarithmic terms 
explicitly displayed in eqs. ( |C.3| ) and ( p.lO| ). 

3.2 Generalized scheme ansatz for Fj^ 

From the results of the last section one can now write an ansatz for F^^, also taking into 
account the most general RS dependence to second order in a straightforward manner. 
After some algebra: 

expl^iill i * ,^e./""' X 



where 



TV 1 

5. = f^^^, (3.6) 
27r^ 7o - bo 



Wan 1 



2bo 16 



|(7o - bo) - 27r2(7; - 61) - 2(70 - 60) In a 



(3.7) 



/3:(a,A7i,5i) 



(2feo)2 



/(2)(A7i) + 270(70 - bo) In^ a + In a 



-^«-(«, A71) + ji^Cii^o - bo) - 47r2(7i - 60) 



(2feo) 



(3.8) 



where f^^^ is a complicated expression given explicitly in Appendix C. 
An important point is that the subtraction function H{g'^) in (|3.3| ) gives no extra contribu- 
tions to the y integral after introducing the contour integration as described in section 

In the general case the renormalization of dim= 4 (gauge-invariant) operators involves a mixing 
matrix for the operators m{qq), m"^ 1, and G^'^G^^. However, the entry for m{qq) only involves m* 1 
and m{qq) itself |Q. As for F^, the only gauge- invariant operator of dim ~2 is m^l. 
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2 and Appendix B, the purely perturbative subtraction gives an analytic function of y. 
Obviously however, the actual values of the perturbative coefficients in ( |3.7| ), ( ^.8| ) are 
affected by the subtraction. As previously mentioned, the unambiguous determination of 
f^^ in ( ^.8| ) involves, in a given (minimal) RS, the knowledge of the three- loop coefficient 
of 1/e, obtained from the perturbative expression eq.( P.3|) of Appendix C. Note also that 
the dependence upon A71 in 



is only through 7f — > 7^^, and the dependence on Bi 
only appears in the second order term, (3^, in contrast with the mass case. As well there 
is no dependence upon B2 or A2 at this (second) order. All these properties are of course 
consequences of the perturbative expression starting at one-loop but zero*'' order. 



The discussion of the previous section on the analyticity domain of ( p.5|) with respect to 
A71 is identical, since the branch cuts are determined by the very same relation defining 
F, eq. ( p.34| ). One can thus proceed to a numerical optimization with respect to the 
RS dependence, along the same line as in section 2.5. The only difference is the Pade 
approximant form to be used: from eq. ( |3.5| ), for m " — > the simple pole behavior is now 
given by the second perturbative term, with coefficient a^, while the first perturbative 
term, 1, gives vanishing contribution for m" — > 0. A Pade approximant taking into 
account these properties, and using the full information as contained in a^, (3'^ is 



'\2 



F' 



{FJ 



1 + A/(F 
1 + p/F' 



which gives the result 



IT, Pade \ 



m 



OA. 
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(3.9) 



(3.10) 



where the dots denote the RS dependence. The region around the optimum of ( 3.1CI| ) 
with respect to the RS parameters is illustrated in Fig. 2. In that case again we found 
a non-trivial flattest extrema minimizing the curvature (plateau), corresponding e.g. for 



the Uf = 2 case to a, 



opt 



4.21, B, 



opt 



—0.017. The optimum values are 



F^^-^-(opt) ^ 0.55 A^(2) (0.59 A^(3) 



(3.11) 



for Uf = 2 (3) respectively. This is also summarized in Table 1, section 5. From (|3.11|) , 
with ~ 92 MeV (and neglecting any explicit symmetry breaking effects due to the 
non-zero u,d,s masses), one thus obtains A;^(2) ^ 167 MeV and A;^(3) ~ 156 MeV. 
Note that these relatively low values are more in agreement with earlier estimates of 
35| , p8| than with the recent experimental measurements of as{Mz) at LEP |5l 



the latter giving a larger Ajj^{3) if evolving from the Mz scale down to very low Q with 
naive perturbation theory. This does not in principle exclude our results, due to the still 
large uncertainties on AjYg{2) and AjYg{3) at very low energies. Moreover, as pointed out 



in ref. |^0|, it could always be that new physics contribution, e.g. supersymmetry, affects 
the running of as between the Mz scale and low energy. 

Finally, like in the mass case, we also give for comparison the (a priori ambiguous) corre- 
sponding value of (|3.10|) in the original MS scheme: 



F^P'^^^MS) ^ 0.40 A^(2) (0.64 A^(3) 



(3.12) 
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4 {qq) 



4.1 Renormalization group invariant m{qq) ansatz 
The (bare) quark-antiquark condensate is most conveniently defined as 

(qq) = lim TrS{x) (4.1) 

where S{x) = i < 0|Tg(0)g(x)|0 > is the quark propagator While the above definition 
( [4.1| ) is gauge-invariant, as is well known (qq) is not separately RG invariant, but m{qq) 
is [Q, which is thus the appropriate quantity to consider when applying the formalism 
developed in section 3.1. Once subtracting the divergences remaining after mass and cou- 
pling renormalization, by following a construction similar to what was done for F^r, we 
obtain a finite ansatz for rh{qq). An awkward situation is that {qq) cannot be directly 
accessed, and has to be extracted from tiny explicit symmetry breaking effects due to 
m 7^ 0. This is of course a well-known problem, not specific to our construction. We 
shall discuss below in section 4.2 a possible way of extracting a {qq){^Ji) value from our 
construction. 

The perturbative expansion up to two-loop order for rh{qq) for m 7^ 0, calculated first 
in ||6^ and independently by us, is given explicitly in Appendix C. Note that the three- 
loop order is not known, which is rather unfortunate since, according to the derivation 
in section 3.1, it implies that one only knows unambiguously the first coefficient of 1/F 
in the m{qq) ansatz. To nevertheless give a more accurate estimate, following the usage 
in such a case, we take into account in our calculations the known RG dependence at 
(9(1/F^) 0. Again, one assumes that the optimization with respect to the RS parameters 
a, A71 etc (whose exact dependence at 0{1/F^) is also known) partly takes into account 
this ignorance on higher order terms. 



The resulting expression, including the full RS dependence, reads 



m{qq) 

MS 



exp[^^J 2bo 




1^ 



r — y ^ 

2iTc J y 



i y/m 



1 



+ ■ 



X 



(4.2) 



with 



Hm) 



Nr. 



47r2 270 - 60 ' 



-(27o-M-vr'(27l-&i)-2( 



^0 



270) In a 



(4.3) 
(4.4) 



"^^ An equivalent definition is from tlie derivative with respect to mp of the vacuum energy. Note that 
in the m ^ hmit, there is an arbitrary phase in the definition (O): the sign of {qq) is thus fixed a 



posteriori, by requiring that for real m > 0, {qq) < 0, to be consistent e.g. with the Gell-Mann Oakes 
Renner ||6l|| familiar relation. In what follows we shall designate by {qq) only the magnitude of the quark 
condensate. 

This allows at least some comparison with the true first order estimates, defined as the expressions 
obtained in the limit foi = 71 = 0. 
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and 

.(2) 



(a, A71, 5i) = [/(S)(A7i) + 470(270 - h) W a - 10^^ In 

- ffa^,-,)(a, A71) + (^Ip [f (270 - h) - 4vr2(27i - 61)] , (4.5) 

where /|^^^ is given in Appendix C. 

4.2 Explicit chiral symmetry-breaking corrections 

Clearly what we are really interested in is the value of the condensate in the exact chiral 
limit, m — >■ 0. Closely related to the latter problem, it is desirable to obtain a connection 



between our construction and the familiar Gell-Mann-Oakes Renner (GOR) relation |61 
e.g. in the SU(2) case, 

- (m„ + md){uu) - Flml + ^(m^ J, (4.6) 

where {uu) = (dd) for exact SU(2) isospin symmetry, ~ 92 MeV and are the pion 
decay constant and mass respectively. 

By definition in ( [4.6| ) rriu^d are the current masses, which break the chiral SU (2) l x SU (2)r 



symmetry explicitly. In contrast, as amply discussed, in (4^) m is an arbitrary param- 
eter, destined to reach the chiral limit m " 0. Accordingly, m ^ for m " 0, so 
that one presumably expects only to recover the result m{qq) for m" 0. This 
is actually the case: although ( |4.2| ) may potentially give a non-trivial result in the chi- 
ral limit, typically the simple pole residue (~ 26o(2C)~'" 6{qq) a(^qq){a), upon neglecting 
higher-order corrections), once we require extrema of this expression with respect to RS 
changes (using for the m'" — > limit a Pade approximant similar to the one for F^), we do 
not find non-zero extrema. Therefore, such a result is not conclusive regarding the actual 
value of (qq) [fl) , although it may be considered a consistency cross-check of our formalism. 

Now, one possible way of extracting (qq) (fi) from the ansatz (|4.2| ) is the following. Let 
us introduce a small explicit chiral symmetry breaking mass term, 

- mo^exprni , (4.7) 

to the basic Lagrangian (|2.3| ). After renormalization, and carefully following the different 
steps as indicated in section 2 and Appendix A, one can show that it amounts simply to 
the following substitution into the integrand of the different previous ansatzs: 

y-^y + ml^p , (4.8) 

except in the factor (dy/y) e^/™", which remains unaffected, and where m'^^p is related to 
rhexp in exactly the same way as m" is related to m, see eq. ( |2.3tj| ). Substituting ( |4.8| ) into 
( ^^ , one can now expand the right-hand side of ( ^^ , firstly for small {y + rn'^^p), and 
take in the resulting expression the true chiral limit, y ^ 0. By finally subtracting rh{qq) 
on both sides it gives an expression which depends on m'^^p only, in the limit m" 0. 
More precisely one obtains after simple algebra 



exp 



(qq) = Urn { (m + mexp){qq) - m{qq)} 
m" 
c-B AAi 



Ai^ (2C)-^^5(,,) 2bo exp[^] {mlpf/^ + Oim^J'/^ , (4.9) 



26, 
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in terms of the quantities already defined after eq. ( ^^ ) , including the RS dependence Q. 



Using the relation between mg^.^ and rfiexp in ( |2.33D , one can now cancel rhe^p on both 
sides of ( f4.9| ) and extract a value of {qq). Note that the dependence upon rhexp = ^expifi) 
does not completely cancel out, since A{A'yi) 7^ 1 in a general RS. Indeed in ([4.9|) there 
is an additional implicit dependence on the scheme, through the relation between m"^p 
and rhexp in (|2.33| ): this is expected, otherwise since m" oc rh/Ajjg one would obtain 



(qq) oc pure number x A|j^, which we know cannot be the case, due to the inherent 
/i-dependence of (qq). Consequently, in contrast with the dynamical mass and F^, it 
would not make much sense to invoke the PMS with respect to the scale dependence a. 



Therefore, we shall fix ( [4 .91) to its MS scheme value (i.e. a = 1, A71 = Bi = 0), which 
indeed makes comparison with other works easier, since the quark condensate is generally 



given in the MS scheme (and typically at p, = 1 GeV) ||35|, Although we do not 
try thus to numerically optimize expression ( |4.9| ), one should remember that the ansatz 
( [4.2|) is assumed to be already optimal with respect to the variational mass parameter, 
m". Indeed, the only remaining contribution from ( [4.2|) in (|4.9| ) at first order in m"xp 
is the simple pole, but with residu given by the RS independent zeroth order terms, as 
a consequence of the factor in ( [4.2| ). Accordingly, potentially non-trivial extrema, 
which may have appeared from the RS dependence at the next orders (as was the case for 
Mq and F^) are washed out. Taking into account the next order m'exp terms in ( [4.9| ) one 
may eventually obtain non-trivial extrema, but those small corrections are not expected 
to drastically change the results given below. 



A more serious trouble when fixing the RS to MS is that the MS scheme value 
A71 = is such that the extra singularities occur at Re[y] > 0, therefore rendering the 
expansion around the origin y ~ a priori completely ambiguous, as explained in section 



2. However, in eq. (|4.9| ), y and m"xp by definition independent parameters, so it 
turns out that the extra singularities, as given by the zeros of dy/dF, are independent of 
m"^p (equivalently, because the extra singularities are due to the Jacobian of the F ^ y 
transformation, and the dye'^^^" term remains unaffected by niexp)- We thus assume that 
the corresponding ambiguities cancel out in the difference in expression ([4.9| ). Although 
we were not able to prove this statement rigorously, at least it is easily checked that the 
rh{qq) part in ([4.9| ), which contains the singularities, do cancel in the first order expansion 



in m''. Within this assumption, we thus obtain the MS result: 



exp 

Ml(/.) ^ 0.647 l^^X'"' ; 0.614 [^^r'\ (4.10) 

where the /i dependence is only due to the remaining dependence on m(/i) = l/2{fhu+m(i) 
or rh{ji) = 1/3 (m^ + fhfi + m,,), for n/ = 2ornj = 3 respectively. Accordingly to extract 
a {qq){fi) value, we obviously have to use as input the explicit quark (current) masses 

^^The (1/^) 7^ 1 power of m^^p cx fh/Kjjg in the second line of ( [4.9[ ) may seem awkward 
at first. It is consistent however with the known perturbative behaviour of (gg)(/i): from ( |4.9| ), 
{qq)[^i) ~ m2''o/T°-i ln(/i/A-pg) for fi > ^.Jlg, so that since m ~ [ln(/i/A;p;g)]"''°/^''° one recovers 
(99)(A') — [ln(/i/A-^)]'''o/^'"' asymptotically, i.e. the condensate runs as inversely as the mass, as it 
should. 
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m(/^). Since the only available values of the light quark masses are extracted from a 
QCD sum rule determination of (ss) and using the mu,d/^s mass ratio (and 



moreover assuming the validity of the GOR relation), it is to be considered more as a 
consistency check than a truly independent prediction. Note however that ( 4.1(j| ) is only 



weakly sensitive to the explicit breaking masses m, due to the small power coefficients 
involved. 

Taking the latest results for the running masses at 1 GeV, as collected in ref. |^ ^ 

m(2) = i(m„ + md) ~ 7.2 ±2.3MeV (4.11) 

and 

m(3) = -(m„ + rrid + m,) ~ 63.1 ± 9.1MeV , (4.12) 

for Uf = 2,3 respectively, this gives 

(gg)i§(lGeV) ~ (0.500 - 0.525) Ajg^(2) ; (0.602 - 0.606) Aj^{3) . (4.13) 

Alternatively, observing that in fact in ([4.10|) the power of rh/Aj^g only depends on the 
RS parameter A71, via the quantity A defined in (|2.35|) , there may be a way of extracting 



{qq){^) in the exact chiral limit, m — 0, by simply choosing a RS such that A = 1 exactly. 
This occurs for A71 = 31/(2887r'^) and A71 = 5/(1927r^) for nj = 2 and 3 respectively, 
and gives 

{qq)f_^,{lGeV) ^ 0.52 A^(2) , 0.58 A^(3) . (4.14) 

Note that the results ( |4.13|) and ( |4.14|) are fairly consistent, given the errors on the quark 
mass rh estimates and the presumed intrinsic error of our method. The not so good 
agreement for n/ = 3 may be attributed to the fact that we have neglected the mass 
of the strange quark in our (chiral symmetric) estimate of Ft,/Kj^ and {qqy^^/Ajjg 
(which is probably not a very good approximation, as indicated e.g. from the fact that 
Fk/Ft, ~ 1.23, experimentally.) 

We thus obtain, within our approach, rather small {qq) values. To make our results consis- 
tent with the sum rules calculations , we would need typically Ajj^ ~ 400 MeV, which 



contradicts the results in section 3 for F^^/ A-j^. (The former value of Ajj^ is however 
more consistent with the experimental results from LEP, as{Mz) = 0.117 ± 0.005). 

There are several ways in which our results may be interpreted. First, it is not excluded 
that our way of performing the small expansion in rriexp in (|4.9| ), although mathematically 
well-defined, may simply not be the physically sensible way of proceeding. Our framework 
amply relies on the continuation to very small values of the arbitrary mass parameter m", 
thus it is may perhaps be only capable of determining the zero order of a chiral symmetry 
breaking expansion in a small mass. This is rather unfortunate for the quark condensate, 
which always appears as suppressed by an explicit mass in any physical relation. Besides, 
we anyhow consider the results in (|4.10|) , (|4.14|) on a less firm basis than the ones for Mq 
and Ft^ in ( p.48| ) and (|3.11|) , due to the lack of PMS optimization (and due to the possible 



remaining ambiguities of (f4.2|) in the MS scheme as well 



^^In ( 4.11 ), ( 4.12 ) the values with errors added linearly were simply obtained from the errors quoted 
in M- 
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As previously we can also choose the RS where Re[ycut] = in (4.2) and (|4.9|), which gives results 
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5 Summary and discussion 

Let us first summarize our main results. 



• We have shown that the variational expansion in arbitrary mass m, as developped 
in the context of the GN model 0], can be formally extended to the QCD case. In par- 
ticular we have seen how to reconcile the variational expansion with renormalization, and 
obtained a non-trivial, finite result in the chiral limit for the mass gap, by resumming the 
variational series expansion in x with an appropriate contour integral. 

• Next, we have exhibited the link between the finite variational ansatz and the RG 
solutions in renormalized form, with the identification of specific fixed point boundary 
conditions. This allows to generalize the construction in a most straightforward way to 
the next RG orders and to more complicated quantities as well, such as the relevant order 
parameters of CSB, {qq) and F^. 

• We have studied in some details how the arbitrary RS dependence influences our 
different ansatzs and, in particular, how the extra singularities appearing in the QCD case 
can be moved away by appropriate RS change at the second perturbative order. The re- 
sulting expressions in a generalized RS have been numerically optimized, within a certain 
approximation, due to the complexity of the complete optimization problem. The optimal 
values obtained for Mg/Ajjg and F.„/ Aj^, summarized in table 1, are quite encouraging. 
Actually one should mention that the (unoptimized) corresponding values in the original 
MS scheme, as given in eqs. (|2.50|) , ( p.l2| ), would equally well be in reasonable agreement 
with other estimates in the literature: given the large present uncertainties on the values 
of Ajyj;5(2), Ajg;5(3) it is difficult to assert which results are closest to the largely unknow 
experimental values. However, the MS values are a priori ambiguous due to the extra 
singularities, as discussed in section 2.4; moreover according to the PMS principle adopted 
here, we consider on general grounds the optimized results as being on a firmer basis. 

• Finally, we have also introduced explicit chiral symmetry breaking mass terms in 
our framework, which plays an especially important role in extracting an estimate of (qq). 
It should also be relevant to a more accurate evaluation of Ft,/ Kjjg and Mq/Aj^ in the 
SU{3) case, where the explicit symmetry breaking due to the strange quark mass is pre- 
sumably not negligible. We plan to study those explicit mass effects in a more quantitative 
way in a future work. 

In table 1, we also give the results that are obtained when restricting the different 
ansatz to the first order, i.e. taking the limit 6i = 0, 71 = and neglecting the 
terms in the different expressions for M^, F^, and m{qq) respectively. (In that case, the 
only arbitrariness is the renormalization scale, and the results were obtained from opti- 
mization with respect to a). As well the "zero order" results, i.e. obtained by taking only 
the pure RG dependence are shown. Strictly speaking, those one- loop (and zero-loop) 



in rough consistency with those in ( 4.13| ), ( [4.14[ ). However, due to the lack of RS invariance of {qq), it is 



then more difficult to compare with the quoted results from other methods |3£ 
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Optimized 
quantity 


zero order 
{pure RG) 
Uf = 2, Uf = 3 


1st order 
(optimized) 

Uf = 2, Uf = 3 


2d order 
(optimized) 

'^Z = 2, = 3 




1.0, 1.0 


1.948, 1.948 

^Opt ^ ^ "! ^ ^ 


2.97, 2.85 

aopt ^ 1.42, 2.84 
5f * ^ .12, .06 




0.355, 0.355 


0.359, 0.360 


0.55, 0.59 

Oopt ^ 4.2, 4.77 

5°^* ~ -.017, - .026 


(m(gg)-P'"*'^)i/4 


0.564, 0.564 


0.669, 0.670 


^0 ~ 




0.65 (.'^ 0.61 (.'^ 


- 


- 


^ — (m 0) 


0.52, 0.58 







Table 1: Optimized results at different orders of the perturbative (non-RG) corrections 
for Uf = 2 and rij = 3 as function of the basic QCD scale Aj^. The corresponding values 
of the RS parameters at optima are indicated, when relevant. For M/Aj^ and Ft^/Ajj^ 
the zero order corresponds to neglecting completely the non-logarithmic corrections, and 
taking the simple pole resulting from (first order) RG dependence. For {qq)^/^ / Ajjg the 
"zero" order corresponds in fact to the (unoptimized) formula ([4.9|) . 



results are again a priori ambiguous, being affected by the presence of extra singularities 
which cannot be moved away. We nevertheless found useful to give them for completeness, 
as they were obtained from naively ignoring the ambiguities and taking the simple pole 
residues in the relevant ansatzs. 



Now another problem is how to estimate the error of the method. As may be expected, 
the rigorous theorems in ref. |^| cannot be applied to the present QCD framework. Ad- 
mittedly, our final numerical results thus require a fair amount of mere "trust" , which is 
the usual problem with the PMS, even indeed when applied to perturbative expressions. 
Note however that all our results are rigorously derived, as far as the pure RG behaviour 
is concerned: it is only the purely perturbative non- RG corrections. Mi ^ Q m (|^), 
which are treated according to the PMS. In other words we assume that it is legitimate to 
treat the remaining perturbative correction with a variant of the PMS. In the GN model 
case [H, the value of the curvature with respect to a at the extremum gave a rather good 
qualitative idea of the error, but there is a priori no more definite relationship between 
the curvature at extrema points and the intrinsic error. The knowledge of the exact GN 
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results allowed a comparison of different Pade approximant forms. As a partial qualitative 
cross-check, we have tried a similar comparison in the present case, and indeed found non- 
trivial optima (when those existed) of the same order of magnitude than ( |2.48| ), ( |3.11| ) 



with alternative Pade constructions similar to those explained in ref. |^ Q We do not 
find sensible to make a systematic study in the QCD case since, as above discussed, it 
would not give much more quantitative idea on the intrinsic error of the method itself. 

Concerning the relatively low {qq) values, letting apart the possibility that our method 
simply fails to extract an accurate estimate in the chiral limit, for the reasons discussed 
above, let us consider taking our {qq) results at face value and discuss their consequences. 
In fact, the possibility of a marginally small (or even zero) quark condensate was recently 
raised in ref. [0, Clearly, {qq) 7^ is a sufficient but not necessary condition for 
spontaneous CSB. For instance, as mentioned in section 3, is an equally well-defined 



order parameter. In ref. |^ the authors moreover emphasize that there are at present no 
clear experimental evidence for a realization of CSB through a large quark condensate. 
As stressed there, the familiar GOR relation is only the first order in fhexp expansion, 
assuming the dominance of the {qq) condensate, precisely. On the other hand, one may 
argue that there are some evidences, in particular from the spectral sum rules |^5[| or 
more recently from lattice simulations [Q, of a low-energy QCD picture with a larger 
{qq). On the theoretical side, giving up the dominance of the quark condensate gives a 
rather complicated framework to describe low-energy QCD: although a generalization of 
chiral perturbation theory is perfectly possible [^, it looses a good part of its predictive 



power. Nevertheless it is certainly not unreasonable to conclude that a small quark con- 
densate is at present not excluded from the data. 

Another possibility would be simply that the dominant contribution to quark condensa- 
tion has a very different origin than the mechanism leading to 7^ and dynamical 
quark masses Mq ^ m^. As mentioned in the introduction other non-perturbative effects 
associated with instantons, which are totally ignored here, are likely to play some role in 



the CSB dynamics. For instance, it was argued long ago |T3[ in some simplified picture of 
the QCD vacuum, that the instanton/anti-instanton interaction can give rise to effective 
(non-local) four-quark interactions of the NJL type 0. 

In any event, the present construction may be considered a definite estimate of the con- 
tributions to order parameters that are associated to the non-instantonic vacuum. 
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•^"In the GN model analysis, we could not use a Pade approximant like in ( ^.47 ), which would give 
inconsistent results due to the fact that h'^^ < 0. 
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Figure Captions. 



Figure 1: dynamical quark mass Mg/Ajjg contribution from expression ( [2.47| ) with nj = 2 
versus the scale parameter a, for different values of Bi near the optimal (plateau) re- 
gion. The exact plateau corresponds to Bi = 0.12. A71 is fixed to 0.00437, such that 
RelVcut] ^ 0. 



Figure 2: dynamical F.,^/Aj^ contribution from expression ( p.lO|) with nj = 2 versus the 
scale parameter a, for different values of Bi near the optimal (plateau) region. The exact 
plateau corresponds to Bi = —0.0173. A71 is fixed to 0.00437, such that Re[ycut] — 0. 
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A Renormalization group material 



We set here some definitions and normalization conventions for tlie renormalization group 
results used at different stages of our construction. 

The homogeneous RG operator (i.e giving zero when applied to a RG-invariant quan- 
tity) is taken as 

d d ^. . d . . d . . ^. 

^rr^dji^^^'^d-g-^-^'^'^d^^ ^^-'^ 

where the RG coefficients are given by 



Pig) = /i^ = -bo -b,g'- b, g' 



(A.2) 



/ \ dffl 2 4 6 

lra[9) = -1-=10 9 +ll9 +12 9 +••• 

m ajj, 



(A.3) 



in terms of the coupling constant g, in a minimal subtraction scheme. In QCD [as = 
g'^/{47r)) the coefficients bi and 7^ are known up to the three-loop order|^B[ and read 
explicitly!^, ^ ^ in the MS'-scheme: 



bo 
bi 

b2 



1 . 2 , 



(167r2)2 

1 .2857 



102 - , 



5033 



325 



'167r2)3^ 2 



18 ^ 54 



(A.4) 



1 

1 404 _ 40 
~ (16^M 

1 ,7494 ,320^,^, 4432, 280 2s . , r^ 

72 = — ( h ( C 3 )nf rif) (A.5) 

(167r2)3^ 3^3^^^ 27 ^ 81 -^^ ^ ' 

where is the number of active quark flavors. 

As is well-known only the first two 6j's and 70 are RS independent. For a gen- 
eral perturbative change of scheme (restricted to second order which is sufficient for our 
purpose), 

9'^9' = 9' {l + Ag' + A,g' + ...), 

m^rri =m{l + + B2g'^ H ) , (A.6) 

the modification of the beta and gamma coefficients in the new (primed) scheme reads: 

60 = ^0 , 

61 = 61 , 

^2 = b 2 - AA + (A2 - Al)bo ; (A.7) 

For a review on RG properties, see the excellent textbook by J. C. Collins ]5^ . 
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and 



7o = 7o , 

I'l = 7i + 260^1 - 70^1 = I1 + A71 , 

7; = 72 + 26i5i + 26o(2fi2-52)-2A7;-7oA2. (A.8) 

Accordingly one can set b'2 = and 72 = by choosing appropriately A2 and B2 in ([A.7|) , 
( |A.8|) , for arbitrary Bi and Ai. These expressions are used in section 2, 3 and 4 to con- 



struct generalized RS dependent ansatzs, at the second perturbative order. 

Next we define the leading log (LL), next-to- leading log (NLL) etc, formal series for 
the quantities we are interested in. Let us start with the mass: from purely dimensional 
considerations, we write for the RG invariant physical (pole) mass, 

= mY^a,{-)f^ (A.9) 

p=o ^^ 

where 



a. 



-)^E«P,^N-)F-^ (A.io) 



H is the MS scale and m = m{n) the Lagrangian mass. The ap^o, (p > 1) coefficients of 
[/n(|)]P define the LL terms, the 0^,1 {p > 2) coefficients of [\n{f )]P-^ define the NLL 
terms, etc. Accordingly the app coefficients are the non-logarithmic perturbative terms, 
at order p. Considering now the equation obtained by applying (|A.1|) on (|A.9|) , as giving 
the /i-dependence of for fixed g"^, one obtains formal series for the LL, NLL, etc |5^ : 
-the LL serie: 

-P ap,o = (70 + 26o(p - l))ap-i,o (p > 1; ao,o = 1) (A. 11) 

-the NLL serie: 

(1 - p) ap,i = (70 + 2bo{p - l))ap_i,i + (71 + 26i(p - 2))ap_2,o + 7o(p - l)ap-i,o, 

{p>2) (A.12) 

and so on. 

The universal (RS independent) LL series can be easily resummed as 

M^^ = m(l + 26o/L)"^. (A.13) 

(where L = ln(^)). In contrast the NLL, NNLL etc series are RS dependent, and not 
straightforward to resum explicitly. We have seen in section 2.3 how to find a resummation 
expression for the NLL series which can be checked to correctly reproduce the terms 
formally defined in (|A.12|) to all orders. 

For of a composite operator of naive mass dimension n and depending only on 
m and g, as considered in section 3 and 4 for and fh < qq > respectively 0, it is 
straightforward to generalize the above expressions in ( |A.11| )- (|A.12| ). One thus obtains 
similar LL and NLL series, that we do not display explicitly. It is easily checked that the 
perturbative expansion of the different ansatzs do reproduce the correct formal series. 
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We only consider however composite operators which can be defined as two point functions. 
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B A contour integral resumming the x dependence 



We review here the formahsm originally introduced in 0, H to obtain non-trivial result 
for the resummation of the perturbative series in x, in the limit x — 1. It leads to the 
dynamical mass ansatz in (|2.26| ), (|2.31|) and its generalization for the composite operators 
in section 3, 4. 

Consider the one-loop RG invariant (bare) expression for the mass Mi as given in 
eq. (U): 

Ml = ^ ^ . (B.l) 

[l-6or(e/2)(47r)^/2^2(Mi)-]-o 

Performing the substitution 

mo ^ mo{l - x); g^x, (B.2) 

provides a new quantity Mi(a;). To pick up the order term in Mi(x) = Z^^o '^9^^ 
(having in mind that we are actually interested in the limit x — > 1), a convenient trick is 
by contour integration: 

rW 1 f . .1 1 



Now performing the sum in ( |B.3D exhibits a (1 — x) ^ factor, cancelling the (1 — x) from 
( [B.2|) . This results in the expression 0: 

mJ"^ ^ — /rfxx"(^+^) mo[/o(x)]~^ , (B.4) 
(x ^ 1) 2«vr J 

where the contour is counterclockwise around the origin, and for convenience we defined 
the (recursive) function 

/o(x) ^1-boX glT[-]m-^{l - x)-' (/o)^-o , (B.5) 

dictated from eq. ( |B.1| ). fo{x) has evidently a power series expansion in x, but, less 
obviously, also admits an expansion in (1 — x), as noted by inverting its defining relation 
( [B.5| ). This implies in particular that x = 1 is an (isolated) pole of Mi. 
Provided that no extra singularities lie in the way, one may distort the integration contour 
in ( |B.4D to go around the cut lying along the real positive axis and starting at x = 1. 
Actually, one can go a step further and reach the q ^ oo limit: after distorsion of the 
contour, only the vicinity of x = 1 survives for g — >■ cx), that one can analyse by changing 
variable to 

l-x=-, (B.6) 
Q 

and rescaling mo by introducing mo = rn'^q, keeping uiq fixed as q goes to infinity. One 
finds in place of ( [B.4| ) 

M, = ^I±e'^ (B.7) 

(g^ oo) 2«7r7 v fQ{v)^''o 



•^^In (B.4) there appeared in fact a factor of 1 — a; from which only the last term contributes to 

the integral due to the analyticity of fo{x) defined in (B.5). 
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where now foiv) = 1—bo gQT[^](mQv) "^(/o)*^^'" • Once performing the renormahzation via 
mo = fhZmi Sq = f^^ZgQ^, Ml in ( |B.7| ) is now finite to all orders: 

Mi = — i dve^^ (B.8) 

where the renormalized function 

/ = l + 26o^?4n[(^)/-(W^^o)]. (B.9) 

We have thus shown how to recover finite quantities with a non-trivial x expansion. 
Eq. ( |B.8| ) however only includes the one-loop RG dependence. To make the connection 
with the pole mass one should include the necessary non-logarithmic perturbative cor- 
rections, already present e.g. at the one loop order in (|2.6[). This can be done without 
affecting the contour integration properties, except that the resulting expression of Mi 
has a more complicated structure around t> ~ 0, which can be however systematically 
expanded around the origin in the way discussed in section 2.2. Generalization of the 
previous construction to the next RG order is possible |^ although here we shall use in 
section 2.3 a more convenient construction directly in terms of renormalized quantities. 



C Perturbative results 

In this section we collect for completeness the known perturbative expansions for the 
three quantities of interest, M(m), F^{rh) and fh{qq){fh) respectively, in the MS scheme. 
These determine, among other things, the non-logarithmic perturbative corrections in our 
different ansatz, as well as serve as a cross-check of the LL and NLL expansion properties. 



C.l Quark mass 



The two-loop pole quark mass was calculated in ref.|^, with exact dependence on the 
(current) quark masses running in the loops. It includes LL and NLL terms plus the 
non-logarithmic perturbative corrections. As we consider only the case with nj equal 
mass quarks, the results of ESI take a simpler, entirely analytical form: 



Mr 



m 



l+7o(^-i^)/ + [7o(| + M^' 



K 



;|^-^7o6o-7i)^ + 7^]/ 



(C.l) 



in terms of the (renormalized) current quark mass m, with L = ln(m//i), and 



K = — ln2H TT"- 

9 18 6 



2 C(3) , 3673 ,7r 



+ 



288 



(— —) 



n 



3) 



V. 



Yoi^^r (C.2) 



for Uf (equal mass) quarks (Cix) is the Riemann zeta-function). We refer to |^ for the 
details of this calculation. 
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C.2 



As discussed in section 3, the perturbative contributions to the pion decay constant F^^ 
may be obtained from the axial-vector-axial-vector two-point correlator evaluated at = 
0, where p is the external momentum. But that is formally the very same quantity which 
appears as the (neutral) part of the p-parameter (at p"^ = 0) in electroweak theory 



up to some trivial overall factors. The two-loop pure QCD corrections with exact quark 
mass dependence were first computed in [Q, and recently the three-loop corrections in 



55 1 . The only caution is to convert some of these results, given in the on-shell scheme, to 



the M5'-schemeP^ More trivially, those expressions were calculated for a top mass rrit, so 
that in the present context we replace rrit — >■ m, where m designates the current (equal) 
mass of the light quarks. One thus has, for D = 4 — e 



-Nr. 



— 2 

m 

'2^ 



L + 



167r2 



3e 3 6 



+ 



9' ^2 



167r2 



^-fn/ 132 -fn/ -910/3 + 24C(3) + 16n 
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where L = \n{m/fi) and 



+ 



+ 



9e 



+ f!'L' + fl'L' + ffL + /f ] } (C.3) 



x30 304 32 
ft = nf 



f 

J 7r 



J TT 

- -8C(3) 



136 32 
149 10 



9 



9 



(C.4) 
(C.5) 
(C.6) 



Ifi 51 98 97 

/f = y^yq~ ^^^^^^ + ^^^^^^ ~ ^^'^ + + y^^^) " Y^^^^ + ^^'^ 

+^%(-^ + ^C(3))-^(2 + 12C(3)) (C.7) 



with 



^4 = 16LZ4(-) + - In^ 2 - -TT^ In^ 2 - —n' 
^2' 3 3 180 



-1.7628.. 



(C.8) 



As explicit, ( P.3|) still contains divergences even after mass and coupling renormalization, 
and should be renormalized by and additional subtraction, according to the discussion in 
section 3. Using eq. (|3.4|) one then obtains after some algebra the perturbative coefficients 
of 1/F and 1/F^ respectively, and (3.,^ given in ( |3.7| )- (|3l^ ), where the two- loop term 
f^^ is given as 



/f (M^) 



11527r4(&o + 7o; 



{[455 - 2Anf + 2Atx\-io + 6o)](7o - &o) 



+ 11527r^(62 - 72) + 9607r^(6o - 7o)7i + 23047r%i(7i - 61) + 36(6o - 7o 



(C.9) 



Results in A/5-schenie are also summarized in [Q. Expression in eq. ( |C.3[ ) have been made consistent 
with our definition of e, which differs from the one in p% pa by a factor of 2. 
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C.3 {qq) 



We have computed independently the perturbative expression of the quark-antiquark 
condensate < qq >, with a non-zero mass, to two-loop order. Our results fully agree with 
the expressions first calculated by the authors of ref. l6^. With D = 4 — e, one obtains 



in the MS'-scheme 

- AT \^ r , " , y ^ " , " , r2 



"0 < 5. >U= iV.^I- - L + - + + - + - -L + _)] (C.IO) 



where the index stands for unrenormalized quantities. Like in the case, the expression 
in ( |C.1CI|) has to be further renormalized by a subtraction, which leads to the perturbative 
coefficients of 1/F and as given in (|4.4| ), (|4.5|), where the two- loop /<gg> term is 



given as 



+ 47r26o7i - 127r%7i - 471 + 2471^?] (C.ll) 



•^^We stress again that f^qqy is noi the complete two-loop perturbative term: the latter would need 
the knowledge of the 3-loop order 1/e coefficient, which we arbitrarily put to zero here. 
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Figure 1: dynamical quark mass Mg/Aj^ contribution from expression ( p. 471 ) with nj = 2 
versus the scale parameter a, for different values of Bi near the optimal (plateau) region. 
The exact plateau corresponds to Bi = 0.12. A71 is fixed to 0.00437, so that Re[ycut] — 0. 
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F_pi/Lambda (n =2) 




Figure 2: dynamical F^^/Aj^ contribution from expression ( p.lO| ) with Uf = 2 versus the 
scale parameter a, for different values of Bi near the optimal (plateau) region. The exact 
plateau corresponds to Bi = —0.0173. A71 is fixed to 0.00437, so that Re[ycut] — 0. 
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